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ABSTRACT 

A  new  technique  has  been  developed  to  study  two  dimensional  heat 
transfer  problems  in  gears.  This  technique  consists  of  transforming  the 
heat  equation  into  a  line  integral  equation  with  the  use  of  Green’s  theo¬ 
rem.  The  equation  is  then  expressed  in  terms  of  eigenfunctions  that  sat¬ 
isfy  the  Helmholtz  equation,  and  their  corresponding  eigenvalues  for  an 
arbitrarily  shaped  region  of  interest.  The  eigenfunction  are  obtained  by 
solving  an  intergral  equation.  Once  the  eigenfunctions  are  found,  the 
temperature  is  expanded  in  terms  of  the  eigenfunctions  with  unknown 
time  dependent  coefficients  that  can  be  solved  by  using  Runge-Kutta 
methods.  The  time  integration  is  extremely  efficient.  Therefore,  any 
changes  in  the  time  dependent  coefficients  or  source  terms  in  the  bound¬ 
ary  conditions  do  not  impose  a  great  computational  burden  on  the  user. 
The  method  is  demonstrated  by  applying  it  to  a  sample  gear  tooth.  Tem¬ 
perature  histories  at  representative  surface  locatons  are  given. 

INTRODUCTION 

The  cooling  of  gears  is  an  important  problem  that  has  been  studied  for 
a  number  of  years.  An  early  model  of  oil  cooling  is  given  in  De Winter 
and  Block  (1972).  El-Bayoumy  et  al.  (1989)  expands  on  the  model  of 
DeWinter  and  Block  (1972)  mainly  by  noting  the  importance  of  the 
Coriolis  force  on  oil  cooling  of  gears.  El-Bayoumy  et  al.  (1989)  also 
develops  a  finite  element  model  of  the  gear  tooth.  We  now  believe  that 
the  finite  element  method  has  speed  and  accuracy  limitations  and  have 
abandoned  this  approach  in  favor  of  the  Green’s  function  method 
described  herein. 

In  the  next  section,  a  new  dynamic,  accurate  and  efficient  solution 
method  for  two  dimensional  heat  transfer  problems  in  gears  is  described. 
The  solution  method  consists  of  transforming  the  heat  equation  into  a 
line  integral  equation  with  the  use  of  Green’s  functions  with  unknown 


time  dependent  coefficients.  The  transformed  integral  equation  is  used 
to  obtain  the  dynamic  equations  for  the  time  dependent  coefficients  for 
each  eigenfunction.  In  order  to  obtain  the  eigenvalues  and  correspond¬ 
ing  eigenfunctions,  the  Helmholtz  equation  for  the  eigenfunctions  is  trans¬ 
formed  into  a  Une  integral  equation  by  the  use  of  the  two-dimensional 
free  space  Green’s  function.  The  integral  equation  is  discretized  into  a 
set  of  homogenous  simultaneous  equations.  The  discretized  version  of 
the  eigenfunction  can  be  obtained  by  solving  a  set  of  homogenous 
simultaneous  equations.  The  obtained  dynamic  equations  can  be  inte¬ 
grated  extremely  efficiently.  Therefore,  any  changes  in  the  boundary  con¬ 
ditions  do  not  impose  a  great  computational  burden.  In  the  third  section, 
the  computational  results  and  a  discussion  is  presented. 

FORMULATION 

In  this  section,  an  accurate  and  efficient  solution  method  for  solving  a 
time  dependent  two  dimensional  heat  problems  in  gears  is  developed. 
The  temperature  field  is  expanded  in  terms  of  the  eigenfunctions  with 
unknown  time  dependent  dynamic  coefficients.  The  dynamic  equation 
for  the  time  dependent  coefficients  of  each  eigenfunction  is  obtained  by 
the  use  of  an  integral  equation. 

Consider  a  gear  tooth  geometry  shown  in  Fig.  1 .  The  transient  heating 
of  the  gear  tooth  is  described  by 

—  =  aV2T  (1) 

at 

where  +  ^^/^y^  and  a  is  the  thermal  diffusivity.  The 

boundary  condition  for  the  left  side  and  top  of  the  gear  is  given  by 
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where  the  equation  represents  the  heat  output  to  cooling  oil,  ^  is  the 
vector  coordinate  on  the  boundary,  K  is  the  thermal  conductivity,  n  is  the 
outward  pointing  normal  unit  vector,  h  is  the  heat  transfer  coefficient 
and  Tg  is  the  temperature  of  the  cooling  oil.  The  boundary  condition  for 
loaded  side  (meshing  surface)  of  the  gear  (right  side  in  Fig.  1)  deter¬ 
mines  the  heat  conduction  into  gear  body  and  is  given  by 

Kf  =  F.(I)  (3, 

where  Fq  is  the  heat  flux  into  the  gear  and  the  remainder  of  the 
boundary  is  described  by  the  vanishing  normal  derivative  of  the  tem¬ 
perature,  3T/9n  =  0. 

Dynamical  Equations  for  Transient  Heat  Flow 

Rather  than  solve  the  physical  equations  directly,  as  in  a  finite  element 
method  or  finite  difference  method,  we  develop  here  a  Green's  function 
method  that  reduces  the  two-dimensional  problem  to  a  one-dimensional 
line  integral  over  the  gear  tooth  boundary.  The  line  integral  equation 
yields  the  dynamical  equations.  This  procedure  yields  a  computational 
advantage  over  the  previous  approaches  where  there  is  no  reduction  in 
the  problem  dimensionally  and  the  time  and  spatial  integrations  are  per¬ 
formed  simultaneously. 


The  heat  equation  is  transformed  into  a  line  integral  equation  by  using 
Green’s  theorem,  (Wyld,  1972) 

(5) 

where  ds  is  the  surface  element  and  dl  is  the  line  element.  Substituting 
Eqs.  (1)  and  (4)  into  the  left  hand  side  of  Eq.  (5)  yields 

where  the  vectors  a  and  c  in  the  integration  limits  are  shown  in  Fig.  1 
(a  driven  gear).  Note  that  a  considerable  simplification  has  been  accom¬ 
plished  for  the  line  integral  on  the  right  hand  side  of  Eq.  (5).  The  only 
contribution  comes  from  the  boundary  conditions  described  by  the  non¬ 
vanishing  normal  derivatives,  Eqs.  (2)  and  (3).  Substitution  of  Eqs.  (2) 
and  (3)  and  expansion  of  the  temperature  field  in  terms  of  the 
eigenfunctions  as 

T{x,t)  =  J^il„(t)v„{x)  (7) 

n 

yield  a  dynamic  equations  for  the  time  dependent  coefficients, 
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Figure  1  .—Simulated  gear  tooth  geometry  and  all  the 
boundary  conditions. 


We  assume,  the  following  eigenfunctions,  and  eigenvalues,  are 
known  and  the  eigenfunctions  satisfy  the  Helmholtz  equation, 


1  2  « 
fin  +tXkn  fin  “ 


^^^mmfirn  ■*"®n  "*"^n 


where 


,  =  Jh(4)Vm(4)Vn(^)dl 


(8) 


(9a) 


B„=-Tjh(4)¥n(4)dl  (9b) 

C„=J_  Fo  (|)<|/n(I)dl  (9c) 

The  line  integrals  are  easily  calculated  as  the  boundary  conditions  are 
changed.  Therefore,  the  calculation  burden  for  changing  the  boundary 
conditions  is  minimal.  These  coupled  first  order  equations  in  Eq.  (8)  are 
integrated  efficiently  with  Runge-Kutta  method. 


(v2+k2)v,/„=0  (4) 

inside  the  gear  tooth  with  vanishing  normal  derivative,  =  0.  (The 

method  for  obtaining  the  eigenvalues  and  eigenfunctions  is  given  in  the 
following  section.)  It  can  be  shown  that  these  eigenfunctions  are  or¬ 
thogonal.  The  eigenfunctions  will  be  normalized  to  unity  by 

||vndxdy  =  l. 


Eigenvalues  and  eigenfunctions 

The  method  for  obtaining  the  eigenvalues  and  eigenfunctions  used  to 
effect  the  simplification  of  the  previous  section  is  given  below.  These 
functions  are  two  dimensional  analogs  of  the  trigonometric  functions 
used  in  Fourier  methods. 

The  eigenvalues,  and  the  corresponding  eigenfunctions,  are 
obtained  through  the  use  of  two  dimensional  free  space  Green’s  func¬ 
tion  (Koshigoe  and  Tubis  1989), 
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G(x|x')  =  -iH^l>(k„|x-x'|)  (10) 


where  is  the  Hankel  function  of  the  first  kind  and  the  Green’s  func¬ 
tion  satisfies  the  inhomogeneous  Helmholtz  equation, 

(v^+k^)G(x|x')  =  -6(x-x')  (11) 

where  6(x)  is  the  delta  function. 

Again  utilizing  Green’s  theorem,  one  obtains 

JJ(g(V'2  +k2)v„  -  +k2)G)ds' 


J 


With  the  use  of  Eqs.  (4)  and  (11)  and  vanishing  normal  derivative  bound¬ 
ary  condition  for  the  eigenfunctions  the  integral  equation  is  simplified 


VnW+JVn(^') 


dl'  =  0 


(13) 


This  integral  equation  yields  the  value  of  at  any  location  inside  the 
gear  tooth  when  the  values  of  on  the  gear  tooth  boundary  are  known. 
Now  let  X  approach  a  point,  ^  on  the  gear  tooth  boundary,  then  Eq.  (13) 
becomes 


|i5»45)+p|yn(5')^®n'-o 


(14) 


in  order  to  specify  the  key  features  of  the  gear  tooth  geometry.  The 
cooardinates  of  these  points  are  given  in  table  1.  The  physical  constants 
used  for  the  calculation  are:  the  gear  thermal  diffusivity,  a  =  0.452  ft^/hr; 
the  thermal  conductivity,  K  =  25  Btu/hr/ft/F;  the  oil  temperature,  T^,  = 
200  F;  the  heat  transfer  coefficient,  h  =  0.34  Btu/sec/ft^/F  follows  from 
DeWinter  and  Block  (1972)  and  El-Baypoumy  et  al.  (1989). 

The  heat  flux,  Fq,  along  the  boundary  from  the  location  ^  ^  to  c 

(shown  in  Fig.  1),  is  given  as  a  function  of  the  distance  measured  from 
the  point  ^  (shown  in  Fig.  2).  This  is  the  heat  generated  in  mesh  for  a  1  - 
inch  wide  gear,  with  a  pitch  radius  of  6  inches,  rotating  a  10,000  rpm  and 
transmitting  500  hp.  We  are  interested  in  gear  steady-state  temperatures 
that  take  hundreds  of  seconds  to  reach.  Hence,  the  detailed  temperature 
changes  that  occur  as  the  gear  goes  in  and  out  of  mesh  cannot  be  re¬ 
solved.  Thus,  both  the  heating  and  the  cooling  have  been  averaged  over  a 
complete  revolution  cycle. 


Table  1 . — Coordinate  of 
Points  in  Figure  1 


Figure  1 
location 

X,  inches 

y,  inches 

a 

-0.197 

0.541 

b 

.073 

.795 

c 

.274 

.337 

d 

.318 

.184 

e 

.371 

.184 

f 

,362 

.000 

The  temperature  calculation  was  performed  using  21  eigenmodes.  (This 
was  found  to  provide  adequate  convergence.)  The  results  presented  herein 
were  generated  on  a  486-PC  with  a  total  running  time,  including  the 
costly  eigenfunction  generation,  of  less  than  one  hour. 


where  is  the  contribution  from  the  singularity  in  the  integrand 
(Burton  and  Miller  1971)  and  is  given  by: 

P^  =  inside  angle  at  the  point  (1 5) 

and  P  designates  the  Cauchy  principal  value  integral.  The  eigenfunction, 
Yjj  is  discretized  in  Eq,  (14)  and  yields  a  set  of  simultaneous  equations. 
The  eigenvalues,  k^^  are  determined  by  setting  the  determinant  of  the 
simultaneous  equations  to  zero.  Once  the  eigenvalues  are  determined, 
the  corresponding  eigenfunctions  can  be  obtained  through  the 
simultaneous  equations.  The  formulation  developed  in  this  section  can 
now  be  applied  to  the  gear  tooth  geometry  (shown  in  Fig.  1)  and  the 
calculation  result  discussed. 


SAMPLE  CALCULATION 

The  technique  developed  in  the  previous  section  is  applied  to  the  gear 
tooth  geometry  shown  in  Fig.  1.  Various  coordinates  are  labeled  in  Fig.  1 


Figure  2.— Time-averaged  heat  input  to  gear, 
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Figure  3. — Gear  surface  temperature  at  selected  locations. 


The  eigenvalues  and  the  corresponding  eigenfunctions  are  obtained 
through  the  use  of  discretized  version  of  Eq.  (14).  The  discretization  is 
carried  out  with  40  line  elements  and  on  each  element,  the  second  order 
approximation  is  used  to  represent  the  geometry  as  well  as  the 
eigenfunctions.  The  temperature  calculations  at  three  different  locations 
on  the  gear  tooth  boundary,  (0.0,0,397)  represented  by  the  solid  line, 
(-0.228,0.066)  by  the  dashed  line,  and  (0.174,0.194)  by  the  dotted  line 
are  shown  in  Fig.  3  as  functions  of  time.  The  curve  labeled  T1  is  centered 
on  the  top  of  the  gear  at  location  (0.0,  0.795).  T2  on  the  cooled  side  of 
the  gear  in  the  region  cooled  by  the  oil  at  (-0.228, 0.464).  T3,  the  highest 
curve,  is  on  the  side  heated  during  mesh,  slightly  above  the  pitch  point  at 
(0. 174, 0.592).  Because  both  the  heating  and  cooling  functions  have  been 
averaged  over  a  complete  gear  revolution,  these  temperatures  can  best  be 
interpreted  as  out-of-mesh  temperatures. 

At  the  present  time  the  calculation  method  is  hard  wired  for  a  single, 
but  representative,  problem.  Furthermore,  we  have  not  yet  included  the 
triangular  portion  of  the  gear  extending  to  the  axis  of  rotation.  Generali¬ 
zation  of  the  method  is  planned  now  that  its  application  to  a  specific 
problem  has  been  demonstrated.  Appendix  A  “Thermal  Analysis  of  Spur 
Gears”  provides  the  geometric  formulation  input. 

CONCLUSION 

A  new  technique  has  been  developed  to  study  two-dimensional  heat¬ 
ing  of  gears.  The  computational  advantage  of  this  technique  over 


previous  approaches  using  the  finite  element  method  or  the  finite  dif¬ 
ference  method  results  from  two  features:  First,  the  problem  is  reduced 
from  two  dimensions  to  one.  Second,  the  time  and  spatial  integrations 
are  separated.  Therefore,  when  compared  with  other  methods,  this  new 
technique  can  provide  substantial  improvements  (one  order  of  magni¬ 
tude)  in  computational  speed.  However,  the  benefit  of  the  dimensional 
reduction  is  manifested  not  only  in  the  computation  speed  but  also  in 
the  ease  of  problem  set-up  since  one  is  required  to  deal  with  the  bound¬ 
ary  not  the  entire  two-dimensional  gear  geometry.  The  other  benefit  of 
this  technique,  based  on  the  separation  of  time  and  spatial  integration, 
is  accuracy.  This  technique  takes  the  full  advantage  of  the  spectral 
method  that  has  exponential  solution  convergence.  This  should  be  com¬ 
pared  to  finite  element  or  difference  methods  where  only  algebraic  con¬ 
vergence  is  possible. 

Since  the  new  technique  is  extremely  efficient,  any  changes  in  the 
time  dependent  coefficients  or  source  terms  in  the  boundary  conditions 
do  not  impose  a  great  computational  burden  on  the  user.  This  result  is 
very  important  when  performing  accurate  Scoring  Analysis  in  gears. 
This  allows  the  bulk  (or  blank)  temperature  to  be  accurately  known. 
Furthermore,  the  gear  out-of-mesh  temperature  is  not  a  constant  along 
the  tooth  profile  at  steady  state  running  conditions,  as  is  often  assumed 
by  gear  engineers.  The  method  is  also  more  adaptable  for  use  in  small, 
lubricated,  concentrated  contacts,  such  as  gears,  since  high  resolution 
can  be  obtained  without  using  large  numbers  of  elements. 

Currently,  we  are  planning  to  extend  this  technique  to  a  multiblock 
application  that  includes  the  remainder  of  the  gear  sector  and  that  further 
optimizes  the  computation  accuracy  and  speed. 
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APPENDIX  A 

THERMAL  ANALYSIS  GEOMETRY  FOR  SPUR  GEARS 

by 

Lee  S.  Akin  *  and  Dennis  P.  Townsend  ** 


ABSTRACT 

The  gear  geometry  needed  to  perform  a  complete  thermal  analysis  for 
a  gear  set  is  very  complex  and  has,  to  the  authors  knowledge,  never  been 
published  before  as  one  set  of  equations  needed  for  the  analysis.  The 
thermal  analysis  of  gears  is  a  very  important  subject  in  that  it  can  be  used 
to  determine  the  scoring  failure  criteria  which  includes  the  blank  (bulk) 
temperature  and  the  flash  temperature  as  used  in  the  Blok  scoring 
formula. 

h  T,-Tf+Tb 

The  geometry  for  this  analysis  must  include  the  involute  geometry, 
the  load  variance  when  the  teeth  mesh  is  in  a  single  or  double  tooth 
contact  zone,  tooth  load  sharing  due  to  varying  deflection,  long  and 
short  addendums,  tip  and/or  root  relief.  Also  included  are  the  affect  of 
heat  partitioning  due  to  varying  tooth  contact  sliding  velocities  over  con¬ 
stantly  changing  Hertzian  contact  band  widths,  gear  set  speed  and  its 
affect  on  the  lubrication  regime  (film  thickness  versus  surface  rough¬ 
ness). 

This  geometry  analysis  will  be  used  as  the  input  parameters  to  com¬ 
plete  the  analytical  computation  of  running  gear  temperatures  using 
Green’s  function. 


INTRODUCTION 

This  is  one  more  in  a  series  of  papers  by  the  authors  on  their  continued 
study  of  the  art  of  the  prediction  of  the  onset  of  scoring  or  scuffing  in 
high  performance  gear  drive  mechanisms  and  our  attempt  to  make  it 
more  scientific.  These  studies  have  examined  a  series  of  disciplines  from 
an  interdisciplinary  lubrication  theory  (ref.  1),  to  a  study  of  the  effect  of 
windage  on  the  lubricant  flow  into  high  speed  gear  teeth  (Ref.  2)  and  a 
model  for  lubricant  flow  in  between  gear  teeth  (Refs.  3,  4,  6,  7,  10 
andll,  13, 16). 

In  addition  these  studies  have  evaluated  the  analytical  and  experimen¬ 
tal  spur  gear  temperature  effects  on  operating  variables  (Ref.  5),  gear 
lubrication  and  cooling  studies  (Refs.  7  to  9,  and  12),  investigations  of 
affect  of  transient  (time  variant)  thermal  and  lubricant  boundary  layers 
(Refs.  14  and  15). 

Most  of  the  above  and  all  of  the  thermal  work  was  done  using  finite 
element  methods  which  produce  large  matrixes  causing  slow  computer 
solutions  to  provide  satisfactory  accuracy.  This  geometry  analysis  is 
developed  for,  "a  computer  program  for  the  computation  of  running  gear 
temperatures  using  Greens  function"  provides  a  new  and  unique  solution 


to  the  gear  temperature  analysis  problem.  The  above  is  accomplished  by 
using  special  integration  techniques  newly  developed  to  accommodate 
the  special  conditions  found  in  lubricated  concentrated  sliding  contacts, 
such  as  found  in  high  performance  aircraft  gear  drives.  This  appendix 
describes  the  intricate  gear  geometry  needed  for  this  analysis,  mostly 
available  in  the  literature  but  certainly  not  all  in  one  place.  This  paper 
hopes  to  fill  that  need. 

FORMULATION  OF  GEOMETRICAL  EQUATIONS 

The  equations  for  the  involute  curve  in  rectangular  coordinates  is  shown 
in  equation  set  2  (see  Fig.  1). 


2.  X  =  (cos0  +  0sin0)  R  or  x  =  R  sin0  see  nomencalture 
y  =  (sin0  +  0cos0)R]3  or  y  =  Rcos0 

Rb 
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The  polar/vectorial  angle  0,  in  Fig.  1  may  be  calculated  from: 


(R^-Rbf 


(R'-Rbf 


=  tan<()-<t)  =  inv(t) 


N  N  +  2 

where:  R  =— &Ro 


so  that  the  pressure  angle  ^  in  Fig.  1  is: 


(r^+R^)  r 

4.  0  =  tan  ^  - - —  =  cos  —and  the  roll  angle  is  £  =  0  +  (|)  = 

Rb  R 

tan(|>.  The  important  radii  of  curvature  equations,  some  at  critical  loca¬ 
tions,  are  shown  below  in  equation  set  5,  and  see  Figs.  1  and  2: 

5.  p  =  Rsin(t)  =  eRj,  and  for  the  mating  gear  =  Rj^sintj)  =  e^Rbrn 
The  lowest  point  of  contact  for  the  mesh  is  p^  calculated  from: 


Pc  =Csin(|)-(Ro+Rb)^ 

Pcm  —  (Rom  Rfam ) 


Rc=  Csin(l)-JR^+Rb  +Rb 


where  C  sin  (b  =  L,  the  line  of  action,  and 


Rlstci 


/  / 
/  / 
/  / 

/  / 

/  / 


\  \  \ 

\  \  \ 

\  \  \ 

\  \  \ 

\  \  \ 


^  \_j,\ 


Pb  =  base  circular  pitch - 

—  Base  circular  pitch  =  Pb- 

Pp^l _ 

- Single  tooth  contact  — 

—  Length  of  line  of  contact  - 

—  Length  of  line  of  action  -- 


Figure  2.— Gear  mesh  line  of  action  geometry. 


^cm  V^om'^^bm  “*"^bm 


is  the  lowest  point  on  the  mating  gear.  The  radii  of  curvature  at  the  low¬ 
est  and  highest  points  of  single  tooth  contact  may  be  calculated  from: 

Pl=(R2+Rgp  Ri=(p?+R^p 
Ph  ””  ~^bm]  'I'Pb 

Rh  =(Ph+Rb) 

271 

where  the  base  circle  pitch  is  Pb  =  —  cos<|)  and  the  circular  pitch 


is  The  distance  of  roll/slide  S  along  the  involute  curve 

may  be  calculated  by  integrating  over  the  roll  angle  8,  from  (vogle. 
Ref.  17),  see  Fig.  1: 


ri?  j  R 
Rbede  =  - 


b  .2  £2  _  i^b  /„2  _2\ 

£,  -1-1^2 -El  j 


and  since,  as  can  be  noted  from  Fig.  1,  p  =  R^e  =  R^  tan  (|)  and  £  = 

2_  2  ^ 
thus  Si_2  =  ~ — ~  over  any  arbitrary  portion  of  the  profile  with  sub- 
2Rb 

2  2 
pQ-Pc 

scripts  1  and  2  and  ^c-o  -  —  over  the  whole  profile  from  the  first 

ZKb 

point  of  mesh  contact  on  the  profile  at  to  the  last  point  at  outside 
diameter  S^.  For  example:  from  R^,  to  Rj  to  R  (at  pitch  point)  to  Rj^  to 
R_  at  the  outside  radius.  The  arc  length  along  the  base  circle  £  can  be 


6 


calculated  from  (were  subscripts  1  and  2  are  at  arbitrary  locations  on  the 
tooth  profile)  as  shown. 


7.  Z,_2  =  Rh(E2  -  e,)  =  Rj,  (tan«t)2  -  tan(t)j)so  that  Xj_2  = 
P2  Pi  1  _ 

P2  Pi  and  thus  the  ratio  between  the  distance  S  and 


R 


the  distance  X  becomes  k  =  a  quantity  useful  in  thermal  calcu- 

2Rb 

lations.  The  transition  time,  as  a  function  of  radii  of  curvature  and  its 

rolling  velocity  between  any  two  points  along  the  tooth,  can  be  expressed 
in  equation  set  8  as: 


8. 


dt  dt  ^  sec 


where: 


a  constant  and  since  p  =  eR^,  then  dE  =  codt  and  —  dt  =  — ,  so 

£  p 


dt 


that 


rP2 


'p,  P  P,  -i(p,-P,) 

P2  Pi  1  P2  ^ 

In  —  fQj.  tjjg  gear  to  rotate  from  £ j 

to  £2*  The  critical  dimensions  along  the  line  of  action  are  shown  in  fig¬ 
ure  2  and  described  below  in  equation  9: 


9.  Z  =  ^Rqi  -  R^i^  +  -^Ro2  -Rb2  -Csintl)  is  the  length  of  the  "line  of 

contact  b  -  f '  as  a  subset  of  the  "line  of  action  a  ~  g." 

Now  we  can  calculate  the  width  of  the  Hertzian  band  of  mutual  con¬ 
tact  at  the  mesh  point  from  equation  10:  (Timoshenko  Ref,  18) 


Therefore  the  rolling  velocity  anywhere  along  the  profile  can  be  cal¬ 


culated  from:  Vi 


=  — -  ft/sec  and  ~  777^  fi/sec  for  the  mat- 
114.59  114.59 


^Pi  "*"^mPim 

ing  gear  so  that  the  sliding  velocity  is  calculated  from:  Vis  ~  ^  — 

ft/sec  where  n  =  speed  of  gear  and  n^^  =  the  speed  of  the  mating  gear. 
Another  value  needed  to  calculate  the  coefficient  of  faction  is  the  total 


velocity  from:  = 


np. 


ft/sec. 


114.59 

We  can  now  calculate  the  coefficient  of  friction  as: 


12.  f  =  0.0127  logio 


'"a.lTxlO^WtSec''^ 
FPcpVsV, 


where  F  =  tooth  face  width, 


=  viscosity  in  centipoise  and  =  tangential  tooth  load  in  lbs.  Thus 


we  find  ourselves  in  a  position  to  calculate  the  instantaneous  heat  flux 
q(p)j  as  a  function  of  the  radius  of  curvature  at  the  instantaneous  posi¬ 
tion  “i”  along  the  line  of  contact  per  equation  Set  13. 


13. 


.V.  X  'WtSec(<!))f(npi-n„Pi„) _ _ 

q(Pi  j  = - - BTU/min 


170330 


10.  B  = 


16W„(Ki+K2)PiP2 

F(Pi+P2) 


n1/2 


where:  - 


_l-vi 


7CEi 


and 


and  since:  p  j  +  P2  =  =  C  sincj) 


V  _  1*“V2  y  y 

K2  -  — - —  and  +  f^2  ~  ^ 

7CE2 

and  =  Wj/cos(l)  =  W^seccj),  see  Fig.  3,  and  Ej  =  E2  =  E 

a/2 


B  =  3.19 


W, 


Fcos(t)Csin<|)E 


{P1P2) 


The  rolling  velocity  for  the  gear  and  it's  mating  gear  V2  are  calcu¬ 
lated  from  equation  set  1 1 : 

11.  Vi  =^^i5PL  ft/sec  and  ^2  ~  ft/sec,  so  that  the  sliding 

360  360 

velocity  %  =Vi“V2  =  ^^(niPi -n2p2)  ft/sec  where:  Vj  =V&V2  = 
Vm. 


This  equation  can  also  be  written  in  a  form  more  useful  using  the 
pinion  speed  only: 

Wt  sec((l))npf  /  ^  _  Ng 

q(pi)  =  — 33Q — (Ppi"~mgpgij  BTU/min,  where  the 

gear  ratio.  At  times  it  is  more  convenient  to  calculate  the  hear  flux  from 
the  radius  vector  at  the  instantaneous  contact  points  where: 


WtSec((t))n  f  f  /  \  f  2 

170330 


-bmg!  Rpi 


Rpo  _J^2 


^pLPC 


Pb 


'■go 

^gLPC 

1/2) 


-R 


gb 


1/2 


BTU/min 


CLOSURE 

This  appendix  develops  the  geometry  analysis  for  the  input  for 
the  computer  solution  of  the  thermal  analysis  of  spur  gears  using 
Green’s  function  to  solve  for  gear  blank  surface  temperatures.  This 
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geometry  analysis  provides  the  gear  tooth  geometry  input,  the 
gear  tooth  sliding  distance  parameters,  the  rolling  and  sliding 
velocity  inputs  and  the  equations  for  the  frictional  heating 
developed  during  the  gear  tooth  meshing,  as  a  function  of  the 
location  on  the  gear  tooth. 

Using  these  inputs  the  program  can  then  determine  the  transient 
and  steady  state  temperatures  of  the  gear  teeth. 
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NOMENCLATURE  FOR  APPENDIX  A 


B 

C 


width  of  the  band  of  mutual  contact  at  the  mesh  point 
center  distance  for  mating  pair 


Rl&Rh 


outside  radius  of  the  gear  and  end  of  the  involute 
curve  and  for  mating  gear 

radius  vector  to  lowest  and  highest  points  of  single 
tooth  contact 


F 

f 


Young's  modulus  of  elasticity 
face  width 

coefficient  of  friction 
gear  ratio  Ng/Np 


V 

V. 


velocity  along  the  curve  (involute) 

the  sliding  velocity  at  an  instantaneous  point  -  V2 
in  the  mesh 

the  total  velocity  at  an  instantaneous  point  in  the  mesh 


N 


n 


P 


number  of  teeth  in  gear  and  number  of  teeth  in 

Vi 

rolling  velocity  of  the  gear  and  its  mating  gear 

the  mating  gear 

w,&w„ 

tangential  and  normal  (perpendicular)  load, 

pinion  speed,  rpm,  n^^^  for  mating  gear  and  Up 

respectively 

for  pinion 

X  &  y 

Cartesian  coordinates  of  the  involute  curve  from  its 

diametral  pitch  of  the  teeth 

origin  at  the  base  circle 

Pb 

qRt 


q(Pi) 


R 


circular  pitch  of  gear  and  its  mating  gear  Z 

heat  flux  calculated  using  the  instantaneous  radius 

vector  At^2 

heat  flux  due  to  sliding  friction  at  instantaneous  point  AN 

of  contact  using  the  radius  of  curvature  as  a 
parameter 

e 

radius  vector  to  the  pitch  point  and  Rj^^  for  the  mate 

at  the  contact  point  0 


length  of  the  line  of  contact  as  a  subset  on  the  line 
of  action 

the  time  it  takes  to  rotate  from  to  E2 

virtual  number  of  teeth  expansion  or  reduction  for 
long  and  short  addendums 

roll  angle  on  tooth 

involute  polar  angle  =  inv(|>  =  tancj)  (rad) 


Rb 

Rc 

Rgb&Rpb 


base  radius  of  the  involute  curve  (its  origin)  and  ^^.p 

for  mating  gear 

P 

radius  vector  to  lowest  point  of  contact  from  center 
of  gear  and  R^j^^  for  mating  gear 

Pc 

base  radius  of  gear  and  pinion  respectively 


^gi 

Rgo*Rpo 


instantaneous  radius  og  gear  and  pinion  respectively 
outside  radius  of  gear  and  pinion  respectively 


Pi&Ph 


oil  viscosity,  in  cp 

radius  of  curvature  from  the  base  circle  and  pj^  for 
the  mating  gear  at  contact  point 

radius  of  curvature  at  lowest  or  initial  point  of 
contact  from  base  circle  and  for  mating  gear  at 
contact  point 

radius  of  curvature  at  lowest  and  highest  points  of 
single  tooth  contact 


^LPC2‘^^LPCi  lowest  point  of  contact  for  gear  and  pinion 
respectively 


pressure  angle  of  mesh 
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